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ABSTRACT 

The  method  of  replication  is  frequently  used  by  simulators  to  estimate 
steady-state  quantities.  In  this  paper,  we  obtain  conditions  under  which  this 
method  yields  asymptotically  valid  confidence  intervals  for  steady-state 
means . 
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sicnii- lcancl  and  explanation 


Consider  a  stochastic  system  which  is  stable  in  the  sense  that  it  obeys 
an  ergodic  theorem.  Simulators  are  frequently  interested  in  estimating  the 
steady-state  ergodic  limit.  For  example,  in  a  queueing  system,  one  often 
wants  to  analyze  the  long-run  number  of  customers  in  queue.  One  popular 
method  for  estimating  steady-state  parameters  involves  generating  independent 
replicates  of  the  stochastic  system;  each  replicate  consists  of  the  time 
evolution  of  the  system  up  to  some  large  (but  finite)  time  horizon.  In  this 
paper,  we  obtain  conditions  under  which  this  easily  implemented  method  yields 
asymptotically  valid  confidence  intervals  for  steady-state  means. 


LIMIT  THEOREMS  FOR  THE  METHOD  OF  REPLICATION 


Pater  W.  Glynn 

1.  Introduction 

Let  Y  «  {Y( t )  i  t  >  0}  be  a  real-valued  stochastic  process  for  which  there  exists  a 
constant  r  such  that 

Y(t)  -->  i 

as  t  *  »,  where  Y(t)  «  t  1  /*  Y(s)ds  and  denotes  weak  convergence.  The  parameter 

r  is  called  the  steady-state  mean  of  Y. 

Simulators  are  frequently  interested  in  estimating  steady-state  means.  One  popular 
procedure  for  accomplishing  this,  called  the  method  of  replication.  Involves  simulating 
independent  copies  of  the  process  Y.  The  independence  of  the  copies  (also  known  as 
replicates)  simplifies  construction  of  confidence  intervals,  since  certain  classical 
statistical  procedures  then  become  applicable;  see  LAW  and  KELTON  (1982),  Section  8.6,  for 
a  discussion  of  this  approach  to  steady-state  simulation.  Our  purpose,  in  this  paper,  is 
to  determine  precise  conditions  under  which  the  method  of  replication  is  asymptotically 
valid. 

One  important  feature  of  this  study  is  that  we  allow  schemes  in  which  the  number  of 

replications  is  allowed  to  converge  to  infinity  with  the  run-length  of  the  simulation. 

This  is  done  for  two  reasons.  Firstly,  it  has  been  shown  by  Schmeiser  [9]  that  in  an 

asymptotic  sense,  confidence  intervals  with  shortest  expected  length  are  obtained  when  one 

bases  the  intervals  on  limiting  normal  distributions  (as  opposed  to  Student  -  t 

statistics).  As  our  paper  indicates  (see  Proposition  3.1  and  Corollary  3.9),  the  limiting 

normal  appears  only  when  the  number  of  replicates  is  allowed  to  diverge  to  infinity. 

Secondly,  infinite  replicate  schemes  allow  the  simulator  to  consistently  estimate  a 
2 

parameter  a  (see  Theorem  3.3),  which  is  itself  of  some  independent  interest.  The 

2  — 
parameter  o  measures  the  asymptotic  variability  of  Y(t)  (see  (2.2)  )i  if  Y(t) 
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nt*asures  the  cost  of  running  a  stochastic  system  at  time  t,  then  r  is  the  long-run 


2  — 

average  cost,  and  o  measures  the  extent  to  which  the  average  cost  Y(t)  may  deviate 

2 

from  r  over  the  interval  [0,t].  Thus,  a  may  itself  be  important  in  determining  the 
suitability  of  a  policy  to  be  evaluated  over  a  planning  horizon  of  t  time  units  duration. 


2 .  A  Central  Limit  Theorem  for  the  Method  of  Replication 

Let  {y  :  i  >  1}  be  a  sequence  of  i.i.d.  replicates  of  the  continuous-time  process 
(rn  order  to  incorporate  discrete-time  stochastic  sequences  U  =  {U^  :  n  >  0}  into  our 
framework,  we  set  Y(t)  =  U[t]/  where  [t]  denotes  the  greatest  integer  less  than  or 
equal  to  t.)  A  replication  procedure  is  a  non-decreasing  function 
m  :  R*^  ♦  R*  (R  =  [  0 , 00 )  )  for  which  there  exists  k  :  R+  +  {1,2,...}  such  that: 

(2.1)  i . )  m(t)  k(t)  <  t 

li.)  m(  t)  k(t)/t  *  1  as  t  “ 

(if  m(t)  >  0,  k(t)  can  be  taken  to  be  [t/m(t)j).  We  interpret  k(t)  as  the  number  of 
replicates  Y i » • • • < ( t ) '  each  simulated  to  time  m(t),  associated  with  a  computer  budget 
of  t  time  units.  Relation  (2.1)i.)  says  that  for  each  t  >  0,  the  replication  procedure 
must  not  use  more  time  than  is  budgeted,  whereas  (2.1) ii . )  requires  that  asymptotically  one 
uses  all  the  time  allotted.  The  estimator  r(t)  available  after  t  time  units  of 
computer  time  is  given  by 

k(t) 

t(t)  “Ml)  j,  Vm(t)) 

ehe,"  V  t)  -  t"1  Jjj  Y,(s)ds. 

As  for  the  process  Y,  we  shall  assume  throughout  this  paper  that; 

!«’,/)  i.)  there  exist  constants  r  and  0  *  0  <  00  such  that 

y?  - 

x(t)  ::  t  i  (Y(t)-r)  ==>  ON(0,1) 

as  t  ♦  ”,  where  N(0,1)  is  a  standard  normal  random  variable  (r.v.) 

2 

ii.)  the  process  (x  (t)  :  t  >  0}  is  uniformly  integrable 
in.)  sup  it1,  >:x(t)  I  :  t  *  0}  <  ». 
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Note  that  (2.2)1. )  implies  that  Y(t)  ==>  r  as  t  ♦  00 ,  so  that  the  parameter  r  is 
necessarily  the  steady-state  mean  of  Y.  Also,  (2.2)11.)  is  equivalent  to  requiring  that 

(2.3)  v(t)  =  EX2(t)  ♦  a2 

as  t  -*■  «*j  see  Theorem  5.4  of  BILLINGSLEY  (1).  Conditions  which  guarantee  (2.2i.)  and 
ii.)  are  available  when  Y  is  a  real-valued  functional  of  a  discrete-time  Markov  chain 
(see  Theorems  1  and  3  in  Section  16  of  CHUNG  [4])  or  when  Y  is  a  real-valued  functional 
of  a  regenerative  process  (see  Lemma  5  and  Corollary  9.1  of  SMITH  [10]).  A  sufficient 
condition  for  (2.2)111 .)  is  provided  by  Theorem  A,  given  in  the  Appendix. 

Our  first  limit  theorem  shows  that  under  a  mild  restriction  on  the  replication 
procedure,  r(t)  is  an  asymptotically  normal  estimator  for  r. 

(2.4)  THEOREM.  Suppose  that  m2(t)/t  +  •  as  t  ♦  ">.  Then, 

V> 

t'2  (r(t)  -  r)  =*=>  ON(0, 1 ) 

as  t  ♦  “. 


covariance  matrix  I.  By  the  continuous  mapping  principle  (Theorem  5.1  of  BILLINGSLEY  [1]) 
applied  to  the  continuous  function  h(x)  «  (x^  +...+  x^J/k  (x  =  (x^ , . . . ,x^ ) ) ,  (2.6)  and 
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(2.7)  show  that 


(2.8) 


V,  -  Vs 

m(t')/2(r(t)  -  r)  ==>  a  k  '2N(0, 


1  ) 


as  n  +  “.  But  by  (2.1)ii.),  k  m!t^)/t^  ♦  1,  so  that  (2.8)  yields  (2.5). 

case  2:  suptkft^)  :  n  >  1}  =  “  and  0  >  0:  choose  a  subsequence  t^  such  that 
k(t^)  *  «•  as  n  +  *>.  Note  that 


(2.9) 

where 


Mf  ) 


(t 


1  /  K  l  h 

*  V2  (rtf  )  -  r)  a  \  " 

n  «  n  *-* 


i=  1 


z,  +  8 

in  n 


zin  -  Yn(Yi )  -  Er(t')) 


Y  =  m<t,>1/2/ik<t*)v<m<t,)>>1/2 
n  n  n  n 

a  =  (t,v(in(t,)])1/2/lk(t')ir(t')),/2 
n  n  n  n  n 

Vs 

8  »  (t' )  2  ( Er  ( t '  )  -  r)  . 


Since  Er(t')  =  EY(m(t'))»  it  follows  that 
n  n 


8  -  ( t ' /m2( t' )  )^2  m(  t'  A2  EX(m(  t ' 


)) 


and  hence  (2.2)iii.)  shows  that  the  assumption  m2(t)/t  ♦  •*  forces  6^  to  tend  to  zero. 

Also,  by  (2.3)  and  (2.1)ii.),  a  *  a  as  n  ♦  •*.  To  treat  the  sum  in  (2.9),  we  view  the 

n 

family  of  r.v.'s  {z  :  1  <  i  5  k( t' ) ,  n  >  1}  as  a  triangular  array.  Note  that  EZ;„  = 
in  n  in 


and 


k(t;> 

l  Ezin 
i-1 


Furthermore,  by  Chebyshev's  inequality. 


EZT 


1 <i<k( t ' ) 
n 


p{z2  >  e}  <  - 

in  e 


k(  t'  )e 

n 


as  n  ♦  °",  so  we  conclude  that  {z^}  is  holospoudic  (see  p.  196-206  of  CHUNG  [5)  for 


results  and  definitions).  To  show  that 


(2.10) 


k  ( t 1  ) 
n 

L 

i“  1 


Zin  =  =  >  N  (  0 ,  1  ) 
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as  n  +  00 


we  need  to  verify  Lindeberg's  condition.  Observe  that  for  n  >  0 


3 .  Confidence  Intervals  for  the  Method  of  Replication. 

Theorem  2.4  does  not  readily  yield  confidence  intervals  for  r,  since  it  involves  the 
unknown  constant  o.  The  type  of  confidence  interval  to  be  used  depends  on  the  behavior 

of  k(t). 


(3.1)  PROPOSITION. 


Suppose 


o2  > 


and  k(t) 


k  >  2.  Let 


1  k<t)  -  2 
r(t)  =  k(t)-1  >  -  r <  t) )  • 


T  hen  , 


(3.2) 
as  t 


k/2  (r(t)  -  r)/r1/2(t)  ==>  tk.1 

where  tk_i  is  a  Student-t  r.v.  with  k  -  1  degrees  of  freedom. 


PROOF.  For  any  sequence  converging  to  infinity,  relation  (2.7)  holds.  Applying  the 

continuous  mapping  principle  to  the  function 

+  1/-_  1  £  _  2  -  Vo 

h ( x )  =  k/2  x  (— — —  l  (x.  -  x.  )  )  2  , 

k  K-i  i  k 

-  2 
where  -  (x1  +. . .+  xk)/k  (note  that  h  is  appropriately  continuous  since  o  >  0),  we 

f  ind  that 

k/2(r(f)  -  r)/r/2(t’)  ==>  tv  .  ,• 

n  n  k  i 

tills  proves  (  3.2  )  . 


r  4  , 

i.I.Jl  THEOREM.  Assume  that  the  process  IX  (t)  :  t  »  0}  is  uniformly  integrable.  If 
kit)  *  “  as  t  ►  “  and  m2(t)/t  ♦  ”  as  t  ♦  then 

mftJHt)  ==>  a2 


as  t  *  °» . 


PKr  h  E  .  Note  that 

’,-4>  m(t,r(t) 


,  k(  t ) 

k(t)  ,L , 
i=  1 

,  k(  t ) 

— 1 ) 

k  ( t ) 

i=  1 


m(  t )  I Y  .  (ml  t )  ) 

l 


m( t ) ( Y . (m( t ) ) 

l 


r(t))2 


-  ro(t)(r(t)  -  r)2 
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?  2  2 

By  Theorem  2.4,  t(r(t)  -  r)*  ==>  o  N(0,1)  as  t  -*■  "  (continuous  mapping  with 


hi x )  =  ) ,  so 


m(t)(r(t)  -  r)2  *  — —  •  t(r(t)  -  r)2  ==>  0 


as  t  »  ",  by  (2.1)ii.)  and  the  fact  that  k(t)  ♦  *.  As  for  the  other  term  in  (3.4),  note 


m(t)E(Y. (m(t) )  -  r)  =  v(t) 


1  ■  -  2 

var(rrET  1  “(t)  (Y1(m(t) )  -  r)  ) 

-  tE|x( t>  |  4  -  v2(t)}  . 


Since  {x  (t)  :  t  >  0}  is  uniformly  integrable,  it  follows  that  EX4(*)  is  a  bounded 
function  (see  (5.1)  of  [1])  so  that  (2.3)  implies  that  the  variance  term  (3.7)  tends  to 
zero  as  t  ♦  °°.  By  Chebyshev's  inequality, 


1  -  2 
P^MtT  £  m(t)(Y1(m(t))  -  r)  -  v(t)|  >  e} 


1  kVt>  _  2 

«  var  (j^(  t)  1  m(t)  (Ya  (m(t)  )  -  r)  )  ♦  0 


as  t  *  ",  by  (3.6)  and  (3.7).  Thus, 


1  k(t)  -  2 
— — -  \  m(t)(Y  (m(t))  -  r)  -  v(t)  “>  0 

k(t)  i=1  i 


as  t  *  “.  Combining  (2.3),  (3.4),  (3.5),  and  (3.8)  yields  the  result. 
The  following  corollary  is  immediate  from  Theorems  2.4  and  3.3.' 

(3.9)  COROLLARY.  Under  the  assumptions  of  Theorem  3.3, 

1/,  V-> 

(3.10)  t'2(r(t)  -  r)/(m(t>r<t))'2““>  N(0,1) 


as  t  ♦  "  provided  o  >0. 

Generation  of  confidence  intervals  for  r,  based  on  the  limit  theorems  (3.2)  and 
(3.10),  is  straightforward.  A  condition  guaranteeing  that  {x4(t)  :  t  >  0}  is  uniformly 
integrable  is  given  by  Theorem  B  in  the  Appendix. 
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Appendix 

Let  Y  be  a  non-delayed  regenerative  process  with  regeneration  times 
0  =  T( 0 )  <  T( 1 )  <«..  .  (See  p.  298-302  of  CINLAR  [6])  for  a  discussion  of  regenerative 
processes.)  Assume  that 

(A1)  E(/p( 1 ’  (1  +  | Y( s)  | )ds ) 8  <  »  . 

Then,  SMITH  [10]  proved  that  Y  satisfies  (1.2)i.)  and  ii.)  and  that  the  steady-state  mean 
has  the  representation 

( A2 )  r=E(/£(1)  Y(s)ds)/ET( 1 )  . 

Set  Tk  «  T(k)  -  T(  k-1 ) ,  N(t)  =  max{k  >  0  :  T(k)  <  t} ,  and 

Zi  =  ^T(i-l)  <Y(S)  "  r)dS  * 

THEOREM  A.  Under  (A1),  ( 2 . 2 ) iii . )  is  satisfied. 


PROOF.  It  is  easy  to  see  that 
(A3) 


N(t)+1 

/  Y(s)ds  -  rt  -  l  ZL  -  / 
k-1  *  C 


T(N(t)  +  1  ) 


( Y  ( s )  -  r )  ds 


From  (A1),  it  is  obvious  that  E(  |  |  +  T  ^ )  <  «•  and  thus  Wald's  identity  implies  that 

N(t)+1 

E (  \  Z.  )  »  E(N(t)  +  1 )  •  E2L  . 

k-1 

(see  p.  137  of  CHUNG  [5].)  By  the  definition  of  Zy,  (A2)  proves  that  EZk  =  0  and  thus 
the  first  term  on  the  right-hand  side  of  (A3)  vanishes.  For  the  second  term,  set 
(A4)  a(t)  =  E(j-T(N(t)  +  1 5 (Y(s)  -  r)ds)  . 

A  simple  renewal  argument  shows  that  a  satisfies  the  renewal  equation 

a(t)  =  b(t)  +  (a  *  F) (t) 


where  *  denotes  convolution  and 

b(t)  -  E{/^(1,(Y(s)  -  r ) ds  ,-  T(  1 )  >  t) 

F(t)  *  P{T( 1 )  <  t}  . 

Since 
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N  { t ) -*- 1 

att)  <  e(  l  /T(k_1 j | V( s)  -  r |ds) 


<  E(N(t)  +  1)  •  E(/^(1)jy(s)  -  r  |  ds )  < 


by  Wald's  equality  and  (A1),  it  is  clear  that  a(*)  is  bounded  over  finite  intervals. 
Hence,  by  Theorem  2.3  of  CINLAR  [6], 


(AS) 

where  U(t) 


aft)  =  (b  *  0) ( t ) 

00 

l  P{T(k)  «  t}  -  EN(t)  +  1.  But 
k-0 

b(t)  <  E{/£( 1 1 ( | V t  s ) |  +  r ) ds  ;  T(1)  >  t} 

<  e{T(1)2  /*(1,(|y(s)|  +  r)ds  i  T( 1 )  >  t} 


<  E{T ( 1 ) 2  /q( 1 > ] Y( s ) | as  +  rT(1)3} 


and  thus 


Mt)  <  c1  min{l,  t  2} 


for  some  c1 . 


So,  there  exists  a  and  c2  such  that 


( A6  ) 


1  a ( t ) |  «  |/[0  ,  b( t-s)U ( ds) | 

<  /(t-«.tj  C2  U(dB>  +  /(0,t-a)  C2  t'’U(dS 

<  c2(U(t)  -  U(t-a) )  +  c 2  U(t)/t  i 


the  first  term  above  is  bounded  in  t  by  relation  (1.7),  p.  360,  of  FELLER  [7],  whereas 
the  second  is  bounded  by  the  elementary  renewal  theorem  (see  Theorem  5.5.2  of  CHUNG  [5]). 
From  (A4)  and  (A5),  it  is  evident  from  (A3)  that 

E(/g  Y(s)ds  -  rt) 

is  a  bounded  function  of  t,  proving  the  theorem. 

THEOREM  B.  Assume  (A1)  holds.  Then, 

EX(t)4  *  304 

,  4  . 

as  t  ♦  and  (X  (t)  :  t  >  0)  ig  uniformly  integrable. 
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PhiJOF .  that 


n 

M(n)  =  l  2k 
k=1 

,5  a  martingale.  For  t  fixed,  set  Mn  =  M((N(t)+1)  n),  where  a  b  i  min(a,b). 

Since  N(t)  +  1  is  a  stopping  time,  it  follows  that  {Mn  :  n  >  l)  is  also  a  martingale, 
i  See  Proposition  5.26  of  BREIHAN  (2]  .) 

Let  Dn  -  =  Z^  I  {n  ( t ) +  1  >  n}.  By  Burkholder's  [3]  square  function 

inequality  for  martingales 

“3  N(t)+1  2  3 

:  A 7 )  max  EM  <  c  e(  jl  0  )  =  c&e(  l  2  ) 

1<k<“  3=1  3  j=1  J 

where  ck  =  18  ( k3//2/(  k-1  ) ) k.  Note  that  since  N(t)  <  “  a.s.,  Mfi  ♦  M(N(t)  +  1)  a.s.  as 
r.  ♦  By  (A7)  and  Patou's  lemma  / 

6  N(t)+1  2  3 

!Ab)  EM(N(t)+1)6  <  lim  EM^  5  cgE(  \  2  ) 

n+“  3=1 

2  2 

Setting  Uj  =  2j  -  £Z  ^ ,  observe  that  Minkowski’s  inequality  yields 
N(t)+1  3  N(t)+1 

( A  9  )  e(  l  zf)  =  E(  I  U.  +  E2^-(N(t)  +  1)  ) 

3=1  3  ]=1  3 


.  N(t)  +  1 

=  V3|  I 

3=1 


0. |3  +  E1/3Z^*E1/3(N(t>+1)3) 

3  1  1  ' 


N(t)+1  32  3i 

<  8  max(  E  |  \  U^|  ,  EZ^EfNI  t)  +  1  )  J  . 


A;pl  y.  nq  the  Burkholder  inequality  a  second  time,  we  find  that 

NItl-i-1  N  ( t )  + 1  3/2  N(t)  +  1  2 

e|  „  U  I3  5  c  E(  \  U*)  <  c Jl  +  e(  l  lT)  )  . 

j=1  3  3=1  3  3=1 

L'-iir.d  5  of  SMITH  [10]  proves  that  under  ( A 1  )  , 

N(t)+1  2 

var(  1  U . )  =  0( t )  ; 

3=1  3 

2 

Wi.  ,  i- quality  therefore  implies  that  the  right-hand  side  of  ( A 1 0 )  is  0(t  ). 

,v.  lemma  a  12,  we  conclude  that  the  bound  in  (A9)  is  0(t3),  so  that  (A8)  proves  that 


PROOF.  It  is  well  known  that  for  any  non-negative  r.v.  N,  EN®  is  bounded  by  a  multiple 
of 

00 

(A13)  l  k1"-1  P(N  >  k}  . 

k=1 

Since  {N(t)  +  1  >  k}  =  { T ( k )  <  t} ,  (A13)  is  bounded  by 


(A14)  1  km_1  +  l  !(”■'  P(Ttk)  <  t  -  kw) 


k<2t/u 

k>2t/W 

(U  »  ET( 1 ) 

,  T(k) 

=  T(k)  -  kp ) ; 

comparison  of  the  first  term  in  (A14)  with  the  integral 

of  g(x)  = 

x®-1 

shows  that  it 

is  of  order  0(tm).  For  the  second  term,  observe  that  If 

k  >  2t/w, 

then 

Wk/2  <  Wk  -  t 

so  that  term  is  dominated  by 

(  A15) 

l  km” 1 • P{T(k)  >  Wk/2}  . 

k>  2t/W 

By  Chebyshev's  inequality  and  (Al),  the  probability  in  ( A 1 5 )  is  bounded  by 
ET(k)S/(wk/2)8,  which  is  0(k_<*).  (since  ET(k)8  *=  0(k4),  as  may  be  verified 
algebraically.)  Thus,  (A15)  is  summable  and  bounded  in  t>  (A14)  then  yields  our  result. 

It  is  worth  noting  that  (AD  is  automatically  satisfied  when  Y  corresponds  to  an 
irreducible  finite  state  Markov  chain,  in  either  discrete  or  continuous  time. 


-11- 


References 


1.  BILLINGSLEY,  P.  Convergence  of  Probability  Measures.  John  Wiley,  New  York  (1968). 

2.  B REIMAN,  L.  Probability.  Addison-Wesley ,  Reading,  Massachusetts  (1968). 

3.  BURKHOLCSR,  D.  Distribution  function  inequalities  for  martingales.  Ann.  Probability 

_1_  (1973),  19-42. 

4.  CHUNG,  K.  L.  Markov  Chains  with  Stationary  Transition  Probabilities.  Springer- 

Verlag,  New  York  (1967). 

5.  CHUNG,  K.  L.  A  Course  in  Probability  Theory.  Academic  Press,  New  York  (1974). 

6.  CINLAR,  E.  Introduction  to  Stochastic  Processes.  Prentice-Hall,  Englewood  Cliffs, 

NJ.  (1975). 

7.  FELLER,  W.  An  Introduction  to  Probability  Theory  and  its  Applications,  Volume  II . 

John  Wiley,  New  York  (1971). 

8.  LAW,  A.  M.  and  KSLTON,  W.  D.  Simulation  Modeling  and  Analysis.  McGraw-Hill,  New 

York  (1982). 

9.  SCHMEISER,  B.  Batch  size  effects  in  the  analysis  of  simulation  output.  Operations 

Research  30  (1982),  556-568. 

10.  SMITH,  w.  L.  Regenerative  stochastic  processes.  Proc.  Roy.  Soc.  London  Ser.  A 
232  (1955),  6-31. 


PWG/jvs 


-12- 


security  classification  of  this  race  cH>«n  o*t*  Knitted) 


REPORT  DOCUMENTATION  PAGE 


report  number 

#2789 


4.  TITLE  (ltd  Submit) 


12.  GOVT  ACCESSION  NO. 


Limit  Theorems  for  the  Method  of  Replication 


READ  INSTRUCTIONS 
DEPORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


S.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Summary  Report  -  no  specific 
reporting  period  . 


t.  PERFORMING  ORG.  REPORT  NUMBER 


7.  authorc*; 


Peter  W.  Glynn 


■  -  CONTRACT  OR  GRANT  NUMBER^*; 

ECS-8404809 

DAAG29-80-C-00  41 


».  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Mathematics  Research  Center,  University  of 
610  Walnut  Street  Wisconsin 

Madison.  Wisconsin  53705 


tl.  CONTROLLING  OFFICE  NAME  AND  AOORESS 


I 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 

Work  Unit  Number  5  - 
Optimization  and  Large 
Scale  Systems 


12.  REPORT  DATE 

February  1985 


U.  NUM8ER  OF  PAGES 

12 


MONITORING  agency  name  A  AOORE SS(I(  dlllotonl  tram  Controlling  OWeo)  15.  SECURITY  CLASS,  (ol  i hit  report) 

UNCLASSIFIED 

IS*.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


See  Item  18  below 


A.  DISTRIBUTION  STATEMENT  (ol  Gila  Rap  on; 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  Iho  obmUoet  entered  In  Block  20,  II  dllloront  from  Report) 


«•.  SUPPLEMENTARY  NOTES 

U.  S.  Army  Research  Office 
P.  O.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


National  Science  Foundation 
Washington,  DC  20550 


20.  ABSTRACT  (Continue  on  rowotom  «/de  ff  n«c moemry  ond  Idmntlly  by  block  numbor) 

The  method  of  replication  is  frequently  used  by  simulators  to  estimate 
steady-state  quantities.  In  this  paper,  we  obtain  conditions  under  which  thi 
method  yields  asymptotically  valid  confidence  intervals  for  steady-state 
means . 


on 

WU  |  jan  73 


COITION  OF  I  NOV  AS  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Ohen  Dele  I 


END 


FILMED 

6-85 


DTIC 


